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Abstract 



o . 

psj , The propagation of the excess of quasiparticles and phonons produced by a 

nuclear recoil inside Sn and Zn superheated superconducting granules will be 
discussed. The decay towards equilibrium of the initial disturbance is assumed 
to be a thermal diffusion process described by a set of coupled heat flow equa- 
tions for the effective quasiparticle and phonon temperatures. The solution is 
carried out analytically for a point source located anywhere inside the super- 
conducting granule with the initial energy distributed in both quasiparticle and 
phonon systems. The calculated time delay between the neutron interaction 
and the nucleation of the phase transition will be compared to the time delay 
distributions obtained by irradiating Zn and Sn SSG detectors with a 70MeV 
neutron beam. 



1 Introduction 



Irradiation tests of SSG detectors with a 70MeV neutron beam have been used to study 
the delay in time between the neutron elastic scattering and the phase transition of 
a single granule inside the detector. Typical delay times of ~100ns and ~500ns have 
been measured in Sn and Zn SSG respectively, as it is discussed elsewhere in these 
proceedings 1 . 

After a neutron elastic scattering, the deposited energy is transferred to electrons 
and to the lattice by the recoiling nucleus. The initial disturbance is localized at the 
interaction point and produces an excess of quasiparticles and phonons in a state very 
far from equilibrium where the kinetic energy of the quasiparticles is larger than kT c . 
After a few ps this hot electron gas exhibits a local equilibrium with the phonons and 
a slower relaxation process dominated by heat diffusion takes place 2 . 

The heat propagation in superconducting spheres is discussed in Ref. 3 for an 
ionizing particle depositing energy only in the quasiparticle system as a primary ex- 
citation. These results can not be applied to the case of nuclear recoils where a 
large fraction of the energy is transferred to the phonons as a primary excitation. 



In the present work we will consider the heat diffusion after a nuclear recoil interac- 
tion located anywhere inside the superconducting granule and with the initial energy 
distributed in both quasiparticle and phonon systems. 

2 Thermal Diffusion 

The decay towards equilibrium of the initial disturbance is described by a set of cou- 
pled heat flow equations 3 for the effective quasiparticle (T e ) and phonon temperatures 
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where D e and D p are the thermal diffusivities, r e and r p are the quasiparticle and 
phonon lifetimes and Tj, the bath temperature. Since SSG detectors are made of 
granules embedded in a dielectric material (plasticine or Al 2 3 powder) the heat 
transfer to the bath is due only to phonons crossing the granule surface with a rate 
1/t s . The initial disturbance is assumed to be a point source in f* at the time t 
described by the delta functions 5 e (f* ,t ) and 5p(f* ,t ) with A e and A p accounting 
for the initial energy share between the two systems. 

Equation [I] can be solved analytically using Fourier transforms and the eigenval- 
ues method 4 with the thermal diffusivities and the carrier lifetimes being constant 
in temperature. Due to the spherical symmetry of the problem, the spatial eigen- 
functions of the solution are the even and odd spherical harmonic functions ^m,n,i(^) 
with the indexes n and m referring to the order of the spherical Bessel functions and 
Legendre polynomes respectively. The boundary conditions on the granule surface 
are used to evaluate the positive zeros of the first derivative of the Bessel functions 
(index I). The final expression for the quasiparticle temperature is: 

T e (r,t) = T b + F e (t) (2) 
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with V the granule volume. A similar expression holds for T p (f, t) replacing F e with 
F p . The time dependent terms F e and F p are linear combinations of the initial quasi- 
particle (Tj) and phonon (T p ) temperatures and exhibit a mixture of relaxation and 
diffusion terms. The terms F e (t) and F p (t) are calculated from F e (t,l) and F p (t,l) 
with / = and do not depend on the spatial coordinates. The derivation of Eq. || is 
discussed in Ref. 5. 

At the beginning of the diffusion process the rise in temperature is localized at 
the interaction point with T e and T p proportional to the partition of the initial energy 
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between the two systems given by the conditions F e (0, 1) = F e (0) = T\ and F p (0, 1) = 
F p (p) = Tp. The terms F e (l,t) and F p (l,t) decay exponentially to zero with time 
and at the end of the diffusion process (t = oo) the equilibrium temperatures are 
T e = T p = T OQ with F e (oo) = F p (oo). The energy conservation inside the granule, in the 
case of no heat transfer to the bath (l/r s = 0), impose the conditions: 
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The initial quasiparticle and phonon temperatures T\ and T % p are evaluated from the 
integral of the specific heat in the superconducting state: 

E r f = V f Te C es dT , E r (l-f) = vf Tp C ph dT (4) 

where / is the fraction of the recoil energy E r transferred to electrons by neutron 
elastic scatterings and can be evaluated from Ref. 6. Typical values of / in Sn and 
Zn absorbers are ~0.22 and ~0.4 for recoil energies of 5keV and lOOkeV respectively. 

The lifetimes ratio and the initial temperatures defined in Eq. |3| and Eq. [|, differ 
from the expression used in previous works 3 where it was assumed that r p / r e = C p h/C es 
and = E r /(VC es ) with the specific heats constant in temperature. This approxi- 
mation holds only at high bath temperatures where the relative temperature rise is 
small. At low temperatures where C es <^C p h, the values of r p evaluated from the ratio 
of the specific heats are an order of magnitude bigger than r e . The carrier lifetime 
ratio defined in Eq. |3] depends on the recoil energy and on the partition of the initial 
energy. 



3 Rise in temperature after a nuclear recoil 

In order to compare the calculations with the neutron irradiation tests, we considered 
Sn and Zn granules with diameters 15/im and 30/zm respectively. The quasiparticle 
diffusivity was extrapolated from reference values 7 using the normal state electrical 
resistivity measured 8 on granules similar to the ones used in the irradiation tests. The 
phonon thermal diffusivity is D p <C D e and can not be extrapolated from measure- 
ments. We performed computations of the quasiparticle and phonon temperatures 
either neglecting D p or with D p = D e /100 without obtaining substantial differences 
in the thermal diffusion process. The lifetimes r e were evaluated from 9 at the tem- 
perature for quasiparticle energies two times the superconducting gap. Typical 
values for r e are ~300ns in Zn and ~lns in Sn. The phonon lifetimes r p were obtained 
from Eq. || The rate 1/t s at which phonons cross the granule surface was neglected 
since t s is of the order of few /is. Nuclear recoils inside the superconducting granule 
were simulated considering a three dimensional grid of 72 equally spaced interaction 
points located in the half sphere with the coordinate <fi in the interval 0— tt. For each 



interaction point, the quasiparticle and phonon temperatures were calculated for a 
point on the granule surface located on the equatorial plane at </> = 7r/2. 

The calculated temperatures are plotted in Fig. 1 for a recoil energy of 5keV 
deposited in a Zn granule at the bath temperature of 50mK. Due to the differences 
between the lifetimes in Zn, r e =300ns and r p =50ns, the initial perturbation is trans- 
ferred to the quasiparticles system in a shorter time scale. At the beginning of the 
diffusion process, the temperature rise depends on the location of the interaction point 
inside the granule and temperatures higher than the final value can be reached 
for interaction points close to the surface. At the end of the diffusion, all the tempe- 
ratures converge to T^. In Sn granules, the time scale for the relaxation process is 
shorter because r e ~r p ~ Ins and the quasiparticle and phonon temperatures exhibit 
the same time dependence. 

4 Time delay 

The nuclear recoil interactions were simulated considering a three dimensional grid of 
122 equally spaced points located in the whole granule volume. For each interaction 
point the temperatures were computed with time steps of 4ns from up to 400ns in 
Sn and 1.5/is in Zn. 

The energy threshold of a single granule is related to the magnetic threshold h — 
l—B a /B s h(Tb) where B s h is the granule superheating field and B a is the strength of the 
applied magnetic field. From the phase diagram it is possible to relate the magnetic 
threshold h to the temperature T* needed to fulfill the condition B a — B sh (T*). 

Previous irradiation experiments on single granules with radioactive sources 10 have 
shown that the phase transition tends to nucleate in a small portion of the granule 
surface (nucleation center) and spread afterwards into the full granule volume. This 
effect was found to be less evident in Zn than in Sn granules. Measurements have 
also shown that there is a dependence of the the phase transition field on the crys- 
tallographic orientation of the granule with respect to the applied magnetic field 11 . 
Since SSG are made of a collection of granules, it is quite difficult to define a criteria 
for the occurance of the phase transition valid for all the granules inside the detector. 

To evaluate the time elapsing from the interaction to the nucleation of the phase 
transition, we monitored the quasiparticle temperature on 5 equidistant points on 
the granule equator and recorded the time for which the condition T e > T* was 
simultaneously satisfied on the 5 points. This is a simplified assumption, because 
the nucleation of the phase transition is assumed to be along the full equatorial plane 
neglecting the possibility of having a more localized nucleation center or a combination 
of the two mechanisms. The measured SSG superheating field distributions 12 were 
used to evaluate the probability distribution of the transition temperatures T*. 

To compare with the measurements, the calculated time distributions were weighted 
with the neutron elastic scattering cross section and with the detection efficiency of 
the neutron counter used in the experiments 12 . Despite the simplified assumption 



used to define the occurance of the phase transition, the calculation well reproduce 
the time scale of the experimental results as shown in Fig. 2. It is important to 
note that the measured distribution have a constant time offset due to the electronic 
delays 1 . 

In the case of Zn granules, where the diffusion time is longer, there is a difference 
in the shape of the calculated distributions for recoil energies in the ranges 5-15keV 
and 40-60keV. At low recoil energies the theoretical distribution is shifted toward 
longer times in agreement with the measurements. In Sn, the quasiparticle relaxation 
times are shorter and both the calculated and measured time delay distributions are 
within 200ns. 

5 Conclusions 

An analytical expression to describe the heat propagation inside superconducting 
granules after nuclear recoils was derived. The calculations differ from previous works 
because the initial energy is shared in both quasiparticle and phonon systems. The 
elapsed time from the nuclear recoil interaction to the nucleation of the phase tran- 
sition was calculated for Sn an Zn granules. The calculations are in good agreement 
with the distributions measured in irradiation tests of SSG with a 70MeV neutron 
beam. 
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Figure Captions 

Fig. 1 Quasiparticle (q.p.) and phonon temperatures versus time for 5keV nuclear 
recoils in Zn. The interactions are distributed in a three dimensional grid of 72 
equally spaced points located inside the half sphere = 0— ir. The temperatures 
are monitored on the granule equator at the point <p = ir/2. 

Fig. 2 Calculated (upper) and measured 1 (lower) distributions of the time delay be- 
tween nuclear recoils and phase transitions in Zn and Sn granules with magnetic 
thresholds h—2% and h—1% respectively. The measured distributions have a 
constant time offset due to the electronic readout. 



